Electronic and spin properties of Rashba billiards 
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Ballistic electrons confined to a billiard and subject to spin-orbit coupling of the Rashba type 
are investigated, using both approximate semiclassical and exact quantum-mechanical methods. 
We focus on the low-energy part of the spectrum that has negative eigenvalues. When the spin 
precession length is smaller than the radius of the billiard, the low-lying energy eigenvalues turn 
out to be well described semiclassically. Corresponding eigenspinors are found to have a finite spin 
polarization in the direction perpendicular to the billiard plane. 
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Spin-dependent phenomena in semiconductor nanos- 
tructures have become the focus of strong interest re- 
cently EI • In nonmagnetic systems, intriguing effects 
can arise from the presence of spin-orbit coupling. Struc- 
tural inversion asymmetry in semiconductor heterostruc- 
tures has been shown |3( to give rise to a spin splitting 
of the same type as was discussed in an early paper by 
Rashba 4] . Its tunability by external gate voltages H |g 
has motivated the theoretical design of a spin-controlled 
field-effect transistor |7|. Novel spin properties arise from 
the interplay between Rashba spin splitting and further 
confinement of two-dimensional electrons in quantum 
wires HUES EH or dots El El El El El E3 Spin- 
orbit coupling has also been shown to affect the statistics 
of energy levels and eigenfunctions as well as current dis- 
tributions [3E1- 

In this work, we study a Rashba billiard, i.e., non- 
interacting ballistic electrons moving in finite 2D regions 
whose dynamics is affected by the Rashba spin-orbit cou- 
pling. The Rashba spin-orbit coupling strength can be 
conveniently measured in terms of a wave-number scale 
2fc so , which corresponds to the Fermi-wave-vector dif- 
ference for the two spin-split subbands. Typical val- 
ues for the spin-orbit -induced spin precession length 
-^so = Tr/k so are of the order of a few hundred nanome- 
ters [l| . The relevant parameter characterizing a Rashba 
billiard of size L is k so L. Taking L = 10/im for a typical 
size of quantum dots, the relevant parameter in Rashba 
billiards can be as large as 70. Furthermore, the tunabil- 
ity of the Rashba spin-orbit coupling strength is a con- 
venient tool to induced changes of the billiard's energy 
spectrum without applying external magnetic fields. 

Below we present interesting features of the energy 
spectrum for Rashba billiards, focusing especially on its 
negative energy eigenvalues. We will show that the den- 
sity of states (DOS) is singular at the bottom of the spec- 
trum. This singular behavior occurs independently of the 
billiard's shape and is most striking if the Rashba param- 
eter is large. We have found that for a circular shape, the 



DOS has additional singularities at negative energies. We 
obtain analytic results for their positions. Their corre- 
sponding eigenspinors have a finite spin projection in the 
direction perpendicular to the billard plane, which is the 
direct result of imposing hard-wall boundary conditions. 

Our central quantity of interest is the Green's function 
for Rashba billiards. Having obtained it, we can derive 
the density of states (DOS) g(E) and the smooth count- 
ing function N(E), i.e., the leading Weyl law [3 EH- 
plus correction terms to it. (A discussion of these con- 
cepts for billiards without spin-orbit coupling can be 
found in Refs. 13 13 Hi 13) We provide ana- 
lytic results for circular Rashba billiards and give the 
first two leading terms for arbitrary shapes based on 
the image method of Berry and Mondragon 27]. The 
latter was developed for neutrino billiards, which have 
two-component wave functions and are rather similar 
to the Rashba billiards discussed here. Comparison of 
our analytical results from the semiclassical treatment to 
the numerically calculated exact energy levels of circular 
Rashba billiards demonstrate perfect agreement between 
the two. An asymptotic expansion of spinor wave func- 
tions for negative-energy eigenstates yields a finite spin 
polarization in the direction perpendicular to the billiard, 
in contrast to the familiar result for a 2D plane. 

In the one-band effective-mass approximation the 
Hamiltonian with Rashba splitting in 2D is given by |3l 
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where <J x ,a y are Pauli matrices. This Hamiltonian gov- 
erns the electron dynamics inside the billiard with Dirich- 
let boundary conditions at the perimeter. (See Ref. (3|.) 
7r/fc so is the spin-precession length, which can be tuned 
independently of the system size HUH- In the absence 
of any lateral confinement, the energy dispersion splits 
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into two branches 128 
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where k = yk% + k y . In the range < k < fc so , one 

branch has negative energies bounded from below by 
— A so = —h 2 k 2 /(2m*). The spin splitting is a conse- 
quence of broken spin-rotational invariance. The spin of 
energy eigenstates, which are labeled by a 2D vector k, 
is polarized perpendicularly to k |28j . Hence, no com- 
mon spin quantization axis for single-electron states can 
be defined in the presence of spin-orbit coupling. At 
a given energy E, two propagating modes exist whose 
wave vectors can be found from the dispersion relation 

©: k ± = k T fcso, where k = ^ + k 2 so . 

The Hamiltonian H commutes with the total angular 
momentum operator J z = —ihd v + § a z , where ip is the 

polar angle. The eigenspinors \Xm ) corresponding to 
the two bands of the Hamiltonian (JTJ are given, in the 
representation of polar coordinates r, ip, by 



,(±) \ = 



± J m (k±r) 
J m+1 (k ± r) 
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Here J m (x) is the Bessel function of integral order m 
and E > — A so . Other independent solutions can be 
obtained when J m (x) are replaced by Y m (x), Hm\x), or 
Hm\x). For Rashba billiards with arbitrary shape, the 
eigenstates can be expanded in the basis J2J using linear 
combinations of both ± spinor states. 

To proceed further, we need the free-space Green's 
function for the Rashba Hamiltonian . Using the fact 

that U 2 
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reads 
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Here E can be a complex number. We note that, for 
negative energies, the retarded Green's function contains 
incoming circular waves besides outgoing waves. 

To satisfy the boundary conditions that G(r,r') van- 
ishes at the boundary, the Green's function is decom- 
posed into two parts: G = G M + Gh, where the homo- 
geneous part Gh fulfills (fi-ij) G H =0. This latter 
Green's function is constructed, in the usual way, from 
the eigenspinors (J3J) as follows 
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The coefficients A m , B m ,C m and D m should be chosen 
such that the total Green's function satisfies Dirichlct 
boundary conditions. In general, one gets an infinite set 
of inhomogeneous linear equations for the coefficients, 
which can be solved only numerically. The circular bil- 
liard is a special case for which the coefficients can be 
given analytically (not presented here) for any m. 



be obtained by g(E) = 
where the trace means the 



From G, the DOS can 
— — lim 1) „>Q+ ImTr G{E + irj) 
limit r — > r', integration of r over the area of the billiard, 
and the trace in spin space. The counting function is de- 
fined by N(E) = J® g(E')dE' and its smooth part N(E) 
requires averaging over a small energy range around E. 



We now consider the circular Rashba billiard of radius 
R. Followi ng th e ideas of the systematic method of Berry 
and Howls [25j, we have calculated the first few leading 
terms of N(E). (Details of the lengthy calculation will 
be published elsewhere.) The result is: 
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where the dimensionless energies e = 2m*ER 2 /h 2 and 
£ so = 2m*A so R 2 /h 2 = k 2 a R 2 have been introduced. 
Here K{x) and E(x) are the complete elliptic integrals of 
the first and second kind, respectively, with the same def- 
initions as in Ref.l3(l The first term for both positive and 



negative energies is the contribution from Goo , while the 
remainder originates from Gh- We note that in a com- 
pletely different context, namely for annular ray-splitting 
billiards, a similar Weyl formula has been calculated[3l| 
involving also elliptic integrals. 
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For arbitrary shapes of Rashba billiards, we can also 
determine the area and the perimeter terms of the 
smooth part of the counting function. The largest con- 
tribution to N(E) comes from TrG^, which is always 
proportional to the area A of the billiard: 



N^E) 



+ A so , for E > 0, 



2irtf (y2^V£ + A so , for - A so <E< 0. 
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It follows directly from Eq. Q that, for negative energies, 
the DOS shows a 1 / \[E + A so singularity at the bottom 
of the spectrum E — > — A so . The perimeter term can be 
derived from the generalization of the image method of 
Ref. Hj| using only the free space Green's function. The 
calculation is very much similar to that applied by Berry 
and Mondragon |27| for neutrino billiards and yields 
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valid for all energies E > — A so . Here C is the perimeter 
length of the billiard. The minus sign is a consequence of 
Dirichlet boundary conditions. For zero spin-orbit cou- 
pling, Ni (E) + N 2 (E) coincides with the previously de- 
rived result for 2D billiards [2l| (apart from a factor 2 due 
to spin) . The area term (JJJ can alternatively be derived 
from the classical phase-space integral in the underlying 
classical approach. However, the classical dynamics of 
electrons in Rashba billiards is described by two Hamil- 
tonians [S2l |. which are reminiscent of the two disper- 
sion branches (0 ■ The constant-energy surfaces in phase 
space arc different for the two Hamiltonians, yielding dif- 
ferent contributions to the classical phase-space integral. 
A simple calculation gives then Eq. (JTJ. 

To get better agreement between the numerically ob- 
tained exact counting function and N(E) one should cal- 
culate further terms besides N±(E) and N 2 {E) for large 
fc so . This motivated us to consider more corrections in 
Eq. |JBJ involving elliptic integrals. For the case of circu- 
lar Rashba billiards, the Schrodinger equation is separa- 
ble in polar coordinates. The resulting radial equation for 
both spinor components leads to the secular equation |3^| 
J m (k + R)J m+ i(k-R) + J m (k^R)J. m+1 (k + R) = 0, where 
m is an integer. This equation is invariant under the 
change m — > — m — 1 (Kramers degeneracy). Formal solu- 
tions of the secular equation having zero wave vector are 
excluded since the corresponding wave functions vanish 
everywhere inside the billiard. We obtain the exact N(e) 
from our solutions of the secular equation for different m. 
In Fig. nji, we plot the difference AN = N(e) — N(e) as 
a function of the dimensionless energy e. The difference 
fluctuates around zero, which shows we did not miss lev- 
els (the mean value of AN is a sensitive test for missing 
levels, see e.g., Ref. 0). Fig.^ shows data for approxi- 
mately 54570 levels. Without correction terms in Eq. © 
with elliptic integrals, AN would increase monotonically 
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FIG. 1: In panel a) the difference AN between the ex- 
act counting function and N(e) from Eq. J^J is plotted for 
^/Eso = k so R = 70. In panel b) the exact counting function 
(solid line) and N(e) (dashed line) are shown. In both panels 
dimensionless energies e = 2m*ER 2 /h 2 are used. 



on average, and would predict rs 70 missing levels in the 
energy range plotted. 

In Fig. ^3, the exact counting function is shown to- 
gether with the Weyl formula © for negative energies 
near the bottom of the spectrum — e so . The overall agree- 
ment is good but the exact N(e) shows an additional 
rounded step structure at certain energies e n . This fea- 
ture shows up only for negative energies, although for 
larger energies this is less pronounced. The step struc- 
ture results in large deviations AN at energies e„ and 
concomitant large peaks in the DOS. To see the reason 
for this behavior, it is useful to plot the energy levels as 
functions of m, as shown in Fig. El The curves in the fig- 
ure start almost horizontally at e n , n = 1, 2, . . . resulting 
in large peaks in the DOS at the same energies. 

Using Hankel's asymptotic expression for Bessel func- 
tions with large argument |35| , we were able to derive the 
energy dispersion accurately in next-to-leading order: 
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valid only for negative energies and n = 
1, 2, . . . , int(2y / e^'/7r). For small m,n the above 
expression agrees excellently with the numerics (e.g., 
eo,i is accurate up to 8 digits for e so = 70). It is 
straightforward to obtain corresponding spinor eigen- 
states and calculate their expectation value for the z 
component of spin. Similar to the case of Rashba-split 
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FIG. 2: The energy levels (in units of 2m* R 2 /h 2 ) of circular 
Rashba billiards as functions of m. For a given m levels are 
ordered and first, second, etc. energies are connected as m 
varies. 



eigenstates in ri ngs |36j, but in contrast to that of 
quantum wires (HtHl, it turns out to be finite. We find 



(2m + 1) cos (2y/eZ) + 2(- !)«+'"' 



(10) 

The Schrodinger equation (including boundary condi- 
tions) for circular Rashba billiards is separable in polar 
coordinates, thus integrable. Hence, all the level statis- 
tics should be Poissonian (see e.g. Ref. 13^1. Indeed, we 
have found that the nearest-neighbor level-spacing dis- 
tribution P(s) is Poissonian (not shown here). For other 
shapes, the spin-orbit coupling destroys the integrabil- 
ity. Random Matrix Theory predicts the statistics to be 
that of the Gaussian Orthogonal Ensemble (GOE) due 
to time reversal symmetry Some other intermediate 
distribution has been found for rectangular shape 
and for small k so , reflecting the fact that the rectangular 
billiard without SO coupling is integrable. 

Finally, a few interesting open theoretical problems are 
listed. The Weyl formula is essential to develop a periodic 
orbit theory for Rashba billiards. (For normal billiards, 
see Brack and Bhaduri's book in Ref. l2ll and a theory 
in case of harmonically confined Rashba systems is given 
111 

Ref.H3) 

The Green's function method presented in 
this work would be a suitable starting point to calculate 
observables such as the magnetization or persistent 
currents j3|| in Rashba billiards. 

In conclusion, we have presented a study of electron 
billiards with spin-dependent dynamics due to Rashba 
spin splitting. Semiclassical results for the spectrum 
agree well with exact quantum calculations. We find in- 
teresting properties of negative-energy states, including 
a finite spin projection in the out-of-plane direction. 
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